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ABSTRACT 

During a magnetically quiet interval (1740-1850 UT), on February 14, 
1967, the magnetic-field intensity and energetic electron fluxes (E ~ 0. 3 - 
2*0 MeV) at ATS 1 exhibited coherent modulations having a frequency of 
33. 8 cph (period « 106. 5 sec) and a duration of approximately 40 oscilla- 
tions. The electron fluxes and the magnetic field oscillated in phase. The 
field perturbation reached 8Y (peak to peak) in the direction of the unper- 
turbed geomagnetic field. The transverse component of the field perturba- 
tion was practically zero. The characteristics of the observed oscillations 
appear compatible with those of a compressional (magnetosonic) excitation 
of the outer magnetosphere. The substantially radial normal mode is 
perhaps driven by a bounce-resonant interaction with the 15-keV protons 
that populate the quiet-day ring current. 

The work of Paulikas and Schulz was conducted under U. S. Air Force 
Space and Missile Systems Organization (SAMSO) Contract No. F04701- 
70-C-0059. 



This communication reports the observation of an unusually periodic 
oscillation of the geomagnetic field intensity at synchronous altitude during 
the mid-morning hours (0740-0850 LT) of a geomagne tically quiet day 
(February 14, 1967; EKp = 9). The peak-to-peak amplitude of the field 
oscillation remained ~ 5Y and the oscillation frequency persisted at 
33.8 cph (cycles per hour) for the duration of the event. Although appar- 
ently a rare occurrence, an oscillation that remains so coherent over an 
interval of — 40 oscillation periods should provide some insight, if properly 
exploited, into the underlying dynamics of the earth’s magnetosphere. 

The event was observed by three instruments aboard the ATS 1 
satellite (150°W longitude; ATS LT = UT - 10 hours), the spin axis of 
which is approximately parallel to the earth’s rotation axis and to 
the local magnetic field. The UCLA magnetometer consists of two ortho- 
gonal fluxgates with the axis of each inclined 45° to the satellite spin axis 
so as to permit a vector measurement of the magnetic field B f Cummings 
and Coleman , 1968a, b]. The Bell Laboratories instrument consists of a 
six-element solid-state detector telescope with coincidence circuitry to 
identify particle species and energy. A 20° half-angle collimator looks 
out perpendicular to the satellite spin axis [ Lanzerotti et ah , 1967]. The 
Aerospace instrument consists of three shielded solid-state omnidirectional 
radiation directors with pulse-height circuitry to identify particle species 
and energy. This instrument provides integral flux measurements above 
four different electron-energy thresholds f Paulikas et al. , 1968]. Bec ause 
of poor statistics, data from the proton channels of both particle instru- 
ments are omitted from the study of this event. 



The integral omnidirectional flux Iq into each of the four Aerospace 
electron channels was sampled twice each 40.965 seconds. The two 
measurements were separated by 5, 12 sec; for certain purposes their 
geometric mean is used as a single data point in this paper. A counting 
rate indicative of the differential undirectional flux into each of the 
seven BTL energy channels was measured three times at approximately 
equal intervals in a basic cycle of 81.930 sec. Data points for each of 
the seven channels, therefore, appear at intervals of 27. 3 10 sec (on the 
average). Thus, the Nyquist frequencies associated with data sampling 
are only slightly larger than the oscillation frequency of the observed 
pulsations; Fourier analysis is required to establish the oscillation fre- 
quency precisely. 

The magnetic-field oscillations and particle-flux modulations 
observed in representative electron channels are shown in Figure I. The 
electron fluxes exhibit a modulation approximately in phase with the 
magnetic -field intensity. The oscillatory curves through the data points 
in Figure 1 serve solely to guide the reader’s eye. 

The magnetic-field data shown in Figure 1 are digitized at 15-sec 
intervals and represented in a field-oriented coordinate system. In this 
system, the unit vector £ points in the direction of the unperturbed 

MM 

_ A A 

magnetic field while p and <p lie in a plane orthogonal to £ and point 

MM Mm. 

approximately in the radial and azimuthal directions, respectively. 

The field components in the p<p£ coordinate system are obtained 
from original data expressed (for example) in the XYZ solar-ecliptic 
coordinate system. A low-pass square filter with a cut-off frequency 
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f = 14. 4 cph (0.004 Hz) has been applied to the original data. The 
magnetic-field vector 13 that survives the low-pass filter defines the 
direction of t . The outward-pointing unit vector \ p is defined to lie in 

/S, ^ 

a geographic meridional plane and to be such that p • f = 0. The third 

AM# 

A ^ A 

unit vector is defined by if = s x P . 

Essentially ail of theT magneti^a^fivity at 3378 cph in Figure 1 
resides in the £ component, i.e. , in the component parallel to the 
unperturbed 13-field. The oscillations are virtually absent from the p 
and if components, which are perpendicular to the unperturbed field. 

Thus, the oscillations are purely compressional. Further, if the oscil- 
lations result from a traveling wave, the direction of propagation is per- 
pendicular to B at the near-equatorial point of observation. The Fourier 
spectra of In B^. and of the logarithmic counting rate for a selected elec- 
tron channel are shown in Figure 2 for an interval approximating 1700-1900 UT. 
The spectral densities of all eleven logarithmic electron fluxes and of 
In distinctly exhibit the spectral peak at 33. 8 cph with a full width of 
1. 4 cph at half maximum. This peak is virtually absent in comparable 
spectra of Bp/B^ and B^/B^ , where B^ denotes the low-pass filtered B^. 

Close inspection of Figure 1 suggests that the particle and field (B^) 
oscillations are approximately in phase. Cross-correlation spectra (not 
shown here) between In and In B^ confirm this visual impression and 
show phase agreement to within 5°. The degree of phase coherence is 
dramatically illustrated in Figure 3, where five minutes of data are 
plotted on an expanded time scale; each data point is located at the center 
of the applicable counting interval. The indicated sinusoids are selected 
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to minimize the mean-square deviation from the logarithmic data obtained 
during the five-minute interval. The computed amplitudes and phases of 
the sinusoids for the four omnidirectional electron channels, relative to 
those of In , are indicated in Table 1. It is noteworthy that the degree 
of modulation in all four of these particle channels exceeds that of . 
Further, the modulation is greater at the higher energies. 

Similarly, it is found that the integrated content, or ’’power”, under 
the spectral peak of In invariably exceeds that of In . (As noted 
above, all these peaks have the same spectral width. ) Table 2, obtained 
by comparing the spectral density of each logarithmic particle flux to that 
of In B ^ at the spectral peak (and taking the square root of the ratio), indi- 
cates the relative logarithmic modulation amplitudes characteristic of the 
two-hour interval as a whole. The values thereby obtained for the four 
omnidirectional channels compare favorably with those given in Table 1 
for the five-minute interval. The energy dependence of the degree of mod- 
ulation can in principle be used to infer the spatial structure of the locally 
observed field oscillation. 

The quantitative procedure that makes this inference possible is 
greatly simplified in the case of equatorially mirroring electrons, for which 
J^/B^. remains fixed along a particle trajectory [ Nakada and Mead , 1965; 
Falthammar , 1966; Lanzerotti et al. , 197o], This follows from Liouviile's 
theorem. While the dynamical trajectory of a particle preserves the first 
invariant p, the radiation detectors are sensitive to fixed energy intervals. 
For equatorially mirroring electrons of kinetic energy E = (Y - Om^c^, it 
is not difficult to show that 
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(8 In J d /3 In B^) E 


|l - [( Y + 1)/2Y] (3 In J A /3 In E) B | (d In B^/d In B^) 
+ (3 In J A /0 In B^) e [l - (d In B^/d In B^)] (1) 


where B^ is the field experienced by the representative particle's guiding 
center. In the event of a radial displacement, for example, the value of 
d In B|j/d In B^. should exceed unity; for an azimuthal displacement it should 
not. The formula stated above is of course valid only in the absence of 
drift-periodic echoes f Lanzerotti et al. , 1967; Paulikas et al. , 1968; 

Brewer et al. , 1969; Lanzerotti et al. , 197 1 ]t which indeed were not present 
on February 14. 

The energy-spectral index (-3 In /3 In E)g can be obtained by standard 
methods; i.e., by using the counting-efficiency functions obtained by instru- 
ment calibration f Lanzerotti et al. , 1967], the observed channel-by-channel 
count rates typical of the time interval in question (filtered to suppress the 
oscillations), and a hypothesized functional form (with adjustable constants) 
for J a (E). On a sufficiently quiet day the "radial-gradient" factor 
(3 In /3 In can be estimated by observing the diurnal variation of 

In J A vs In , assuming that the magnetospheric configuration remains 
substantially constant in time. 

Accumulative uncertainties associated with the various observational 
inputs to Equation (1) have prevented a quantitative understanding of Table 2. 
Certain qualitative inferences can be drawn from the analyses, however. 

The energy spectrum during and immediately preceding the micropulsation 
event has been found to resemble those previously obtained at the synchron- 
ous orbit f Lanzerotti et al. , 1967] , i. e. , the spectral index (-9 in J, /9 In E)^ 

i. i 13 

at the higher energies is ~ 5, and therefore is considerably larger in mag- 
nitude than the value ~ 1 found at the lower energies. In other words, the 
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spectrum is quite steep at the higher energies (E > 1 MeV) and quite flat 
at the lower energies (E £ 500 keV). This property would tend to make 
the modulation more intense in the higher-energy channels, as observed. 
However, several values of (9 In /d In in Table 2 would require 

spectral indices ~ 10 in the event that d In B^/d in = 1. Since this 

requirement is not satisfied by the observed energy spectrum, it must be 
concluded that the February 14 micropulsation is not a wave that propa- 
gates purely in the direction of azimuthal drift. 

Thus, it seems appropriate to consider a likely alternative, i. e. , a 
substantially radial oscillation of the magnetosphere in a compressional 
mode. Frequency selection might then reside in the resonant properties 
of the magnetosphere. Cummings et al. [l969] have appealed to the idea 
of Alfvenic eigenmodes of the magnetosphere to account for spectrally pure 
transverse (noncompressional) oscillations having periods ~ 100 sec. 

Those oscillations were also observed at synchronous altitude during 
quiet times, but (being noncompressional) exhibited magnetic activity in 
the p<p plane rather than in the £ direction. The present observations 
suggest that perhaps in this case the magnetosphere (or that portion of it 
that lies beyond the plasmapause) has acted as a magnetosonic resonant 
cavity. 

Except for the implications of Equation (1), the observations would 
have been compatible with the occurrence of a compressional drift wave 
carried by 15-keV protons, which have a total velocity v « 1700 km/sec. 

A wavelength of 2 tt gyroradii and a phase velocity comparable to the 
azimuthal drift velocity of these protons at L = 6. 6 would have produced a 


- 7 - 



wave frequency w/Ztt ~ 3v/4irJLa ** 10 Hz, as observed, where a is 

3 

the radius of the earth. The estimated wavelength ( ~ 1 0 km) would have 
substantially exceeded the gyroradius of a 1-MeV electron (~40 km), 
thereby accounting for the apparently adiabatic response of relativistic 
electrons to the field oscillations. 


In view of the greater-than-expected degree of modulation in sev- 
eral of the electron channels, however, the drift-wave interpretation 
must be rejected in favor of a substantially radial magnetosonic oscilla- 
tion, for which d In B^/d In B^. in Equation (1) can significantly exceed 
unity. In this case the observed frequency can be accounted for by 

postulating a wavelength of roughly twenty earth radii, an ambient field 

3 

~100Y, and a plasma density ~3 protons /cm . The gradual growth of 
wave amplitude with time (Figure 1) suggests a convective plasma 

instability within the magnetosphere rather than an externally applied 

> 

impulse. In this context it is interesting that the bounce frequency of a 
15-keV proton is comparable in magnitude to the oscillation frequency 
observed. Bounce-resonant (15-keV) protons should be present in the 
magnetosphere as a major constituent of the quiet-day ring current. 
Either an inward radial gradient in their spatial distribution or an off- 
equatorial peak in their mirror-point distribution at constant \x (first 
invariant) might be an adequate source of free energy to initiate the field 
oscillation. Unfortunately the available instrumentation on ATS 1 did not 
provide for observation of protons in this relevant energy range. 

Evidence for the magnetospheric scale of the February 14 event is 
provided by the magnetic-field data recorded during the same time period 
by the ground-based station at College, Alaska. College lies close to the 
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foot Of the field line passing through ATS 1. The rapid- run magnetogram 
shows a perceptible oscillation having a period ~ 2 min throughout the time 
interval 1740-1900 UT. Events having this character, however, are quite 
commonly observed at College, and coincidence with a compressional wave 
at synchronous altitude is apparently quite rare. 

Indeed, the present observation at synchronous altitude appears 
unique among compre s sional waves in that it occurred during a magnetically 
quiet period. The three hours preceding the event (1500-1800 UT) exhibited 
Kp = 0, and the three hours following onset (1800-2100) showed Kp = 1 + 
r Lincoln , 1967]. Hourly Dst values during this period were +2, +1, +0, 

+3, -0, and -2 gammas, respectively f Sugiura and Cain , 1969]. 

The oscillations observed at this time were quite different from any 
magnetic oscillations previously reported at synchronous altitude; e.g. , 
storm-associated Pc-5 wave events f Barfield and Coleman , 1970], 
substorm-associated irregular pulsations f McPherron and Coleman , 1970], 
quiet-time transverse magnetic oscillations f Cummings et al. , 1969], 
and band-limited pulsations f McPherron and Coleman , 197 l] . Of the pre- 
viously reported oscillations, only the transverse oscillations described by 
Cummings et al. [1969] occur during magnetically quiet periods. The 
storm-associated Pc-5 wave events are distinct from the event of February 
14, 1967, in that they have a large transverse component. The substorm- 
associated irregular pulsations and band-limited pulsations mentioned 
above are quite irregular and occur only during magnetospheric substorms. 
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Moreover, the present radial oscillation differs in almost every 
respect from the azimuthally propagating drift wave detected by Brown 
et al. [1968] at L « 5 on April 18, 1965, and attributed to the drift 
mirror instability [ Hasegawa , 1969? Lanzerotti et al. , 1969]. That wave 
occurred during a geomagnetic storm, exhibited a much larger oscilla- 
tion amplitude (~40Y peak-to-peak) than the present wave, and broke up 
into turbulence after only a few oscillations. Relativistic electrons were 
apparently heated non-adiabatically in the process. Further, the onset 
of the April 18 wave was preceded by a pronounced diamagnetic decrease 
of the field strength, clearly indicating that p (defined as the ratio of 
plasma pressure to magnetic pressure) was then of order unity. Accord- 
ingly, it was possible irTthat case to invoke a fluid-like instability 
(limited by gyroradius effects) rather than one involving 15-keV protons 
acting as resonant particles. Although observations by F rank [1967] 
suggest that p ~ 1 at the center of the ring current even during quiet 
periods, it would be unrealistic to invoke a fluid-like instability to explain 
the present observation because no apparent diamagnetic impulse precedes 
the onset of oscillation (Figure 1). A low- p condition among particles 
distributed unfavorably for wave growth could perhaps account for the 
unusual duration of phase coherence in the February 14 event. 
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Table 1. Amplitudes and phases of log I«(E) relative to log Bi, 
1759-1804 UT, as illustrated in Figure 3. 


Omnidirectional Electron 
Channel Designation 

Relative 

Amplitude 

Relative 

Phase 

X 1 S 

1 0 (E > 

0.30 MeV) 

2. 11 

+ 1.9° 

1 z ~ 

I 0 (E > 

0.45 MeV) 

4.34 

-6.7° 

h S 

I 0 (E > 

1.05 MeV) 

6.21 

! 

i 

+ 

to 

• i 

o 

o 

la = 

I„(E > 

1.90 MeV) 

6.44 

+ 2.2° 
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Table 2. Observed ratio of logarithmic modulation amplitudes, electron 
flux vs magnetic field, as deduced from Fourier spectra (1700-1900 UT). 


Flux 

Type 

Energy 

(MeV) 

Relative Amplitude 
(8 In Flux/a In 

J l 

0.4 

4.74 

J 1 

0. 6 

6. 45 

J i 

0. 8 

8. 14 

J 1 

1. 1 

8. 74 

J 1 

1.3 

8. 12 

J 1 

1. 5 

6. 7 1 

J i 

1.9 

6. 49 

V 

0.30 

2. 53 

*0 

0.45 

4. 85 

«0 

1. 05 

6. 55 

*0 

1. 90 

5. 62 
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FIGURE CAPTIONS 


Fig. 1: 


Fig. 2 


Fig. 3: 


Observations of electron-flux and magnetic -field oscillations 
at synchronous altitude. Absence of magnetic activity at 
33.8 ±0.7 cph in the p<p plane identifies the oscillation as 
compres sional . Logarithms (LOG) refer to base 10. 

Fourier spectra of representative logarithmic particle and 
field measurements, approximately 1700-1900 UT. Log- 
arithms (LN) refer to base e. For computational convenience, 
the analyzed time interval has a length of 7209. 840 sec (= 264 X 
27.310 sec). The computed point nearest the spectral peak 
therefore corresponds to a frequency of 33.9536 cph ( = 68 X 
0.4993176 cph). The subsidiary peak at 22 cph arises from 
an instrumental effect contained in the particle data. Spectral 
behavior at f < f = 14.4 cph is influenced by application of a 

A 

numerical filter to the field data, as required for defining £ . 

MM 

Five-minute detail of logarithmic particle (Iq) and field (B^) 
measurements, indicating phase coherence and relative 
amplitudes of best-fitting sinusoids. Logarithms (LOG) refer 
to base 10, Scale for log B is expanded by a factor of two 
(relative to the logarithmic electron-flux scales) for clarity. 
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